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Abstract 

In a previous paper (Gayon & Bois 2008a), we have shown the general effi- 
ciency of retrograde resonances for stabilizing compact planetary systems. Such 
retrograde resonances can be found when two-planets of a three-body planetary 
system are both in mean motion resonance and revolve in opposite directions. For 
a particular two-planet system, we have also obtained a new orbital fit involving 
such a counter-revolving configuration and consistent with the observational data. 

In the present paper, we analytically investigate the three-body problem in this 
particular case of retrograde resonances. We therefore define a new set of canon- 
ical variables allowing to express correctly the resonance angles and obtain the 
Hamiltonian of a system harboring planets revolving in opposite directions. The 
acquiring of an analytical "rail" may notably contribute to a deeper understand- 
ing of our numerical investigation and provides the major structures related to the 
stability properties. A comparison between our analytical and numerical results is 
also carried out. 
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1 Introduction 

The simplicity and elusiveness of the full three-body problem have captived math- 
ematicians for centuries. Still today, because of its non-integrability, the three- 
body problem remains tricky. It has been solved for particular cases, by reducing 
the number of degrees of freedom. One of the first major results is the deter- 
mination of families of periodic orbits in the three-body problem under certain 
assumptions (see for instance Henon 1976 or Hadjidemetriou 1976). 

Since the first observations of extrasolar systems, the interest of the general 
three-body problem has been boosted, mainly under the coplanar assumption. 
Moreover, the occurrence of bodies in Mean Motion Resonances (MMR) allowed 
to study the stability of several two-planet systems (see e.g. Hadjidemetriou 2002, 
2008 or Callegari et al. 2004). The secular behavior of two-planet systems has 
been investigated as well (e.g. Henrard & Libert 2008). Because most presently 
known extrasolar planets move on eccentric orbits, Beauge & Michtchenko (2003) 
investigated the eccentric, coplanar three-body problem to characterize the dy- 
namics of such two-planet systems. 

Faced to the detection of multi-planetary systems (one star and 2 planets or 
more), and to the increasing number of degrees of freedom, numerical methods of 
global analysis have been performed in order to study their dynamical stability. By 
combining A'^-body integrators with numerical tools for detection of chaos (notably 
fast indicators or frequency analyses), the exploration of the A^-body planetary 
problem (with > 3) is accessible from numerical methods. As a consequence, 
mechanisms generating stability of multi-planetary systems can be identified and 
explained. For instance, mechanisms involving a MMR as well as an Apsidal Syn- 
chronous PrecessiorQ (ASP hereafter) have been intensively studied (e.g. Lee & 
Peale 2002, 2003; Bois et al. 2003, Ji et al. 2003). A solution involving both 
MMR and ASP describes well the stability of eccentric, compact multi-planetary 
systems, but may not however be unique. We note, for example, that other multi- 
planetary systems have been found to be mainly controlled by secular dynamics 
(e.g. Michtchenko et al. 2006 or Ji et al. 2007). In Gayon & Bois (2008a), by using 
a numerical method of global analysis^ we found a novel stabilizing mechanism 
involving a retrograde MMR in a two-planet system, which means that the two 
planets are both in MMR and revolve in opposite direction^ Let us note that the 
difference between a prograde MMR and a retrograde MMR does not only lie 
in the retrograde motion of one of the two planets. The mechanisms of stability 
related to a retrograde MMR and their underlying properties differ from the pro- 
grade case. Such stabilities are generally more robust (see Gayon Sl Bois 2008a, 
2008b). 

The assumption that two giant planets are in a MMR and revolving in oppo- 
site directions around their hosting star is apparently contradicting to the most 
accepted formation theory of planetary systems, notably to the formation and evo- 
lution of the resonant planetary systems (core accretion mechanism combined by 



' In case of ASP, the apsidal lines on average precess at the same rate. 

^ We used the MIPS (Megno Indicator for Planetary Systems) method based on the MEGNO 
technique (Mean Exponential Growth factor of Nearby Orbits; see Cincotta & Simo 2000). 

^ For instance, the inner planet has a prograde direction while the outer planet a retrograde 
one. In the following, we will also use the terminology of "counter-revolving configuration". 
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a planetary migration scenario). However, in Gayon & Bois (2008a), we present 
two feasible processes leading to planets revolving in opposite directions. The 
first scenario has been introduced by Nagasawa et al. (2008) and consists in the 
combination of a planet-planet scattering, a tidal circularization and Kozai mech- 
anisms. Starting from a hierarchical 3-planet system and considering a migration 
mechanism including a process of planet-planet scattering and tidal circulariza- 
tion, the authors show that close-in planets may be formed. In a few cases, due 
to the Kozai mechanim, one planet enters a retrograde motion. The second feasi- 
ble process imagined by Varvoglis (private discussions) and ourselves is related to 
the capture of free-floating planets. By integrating the trajectories of planet-sized 
bodies that encounter a coplanar two-body system (a Sun-like star and a Jupiter 
mass), Varvoglis finds that the probability of capture is significant, and, morever, 
that the percentage of free-floating planets captured is higher for retrograde mo- 
tions than for prograde motions. As a consequence, it would be possible to find 
one day some planetary systems in counter-revolving configuration. 

In the case of a two-planet system, analytical and numerical methods are com- 
plementary. Besides, the acquiring of an analytical "rail" contributes to a deeper 
understanding of the numerical exploration. As a result, the analytical approach 
permits to study thoroughly the numerical results. By plotting surfaces of sec- 
tions, it provides the major structures related to the stability properties. As a con- 
sequence, in the present paper, we propose a Hamiltonian expansion based on this 
novel and feasible stabilizing mechanism introduced in Gayon & Bois (2008a 
and 2008b). We therefore investigate the three-body problem in the particular 
case of retrograde resonances, adapting the Hamiltonian approach of Beauge & 
Michtchenko (2003) firstly expanded for eccentric, coplanar, and prograde orbits. 
From this analytical expansion, we can correctly derive the expression of reso- 
nance angles in the case of a retrograde MMR. A comparison with the numerical 
results obtained in Gayon & Bois (2008a) is also carried out. 

In the present paper, we firstly introduce the general expressions of the Hamil- 
tonian as well as the canonical variables in the case of two massive planets mov- 
ing around their host star (Section 2). We note that a dynamical system composed 
of two satellites orbiting a planet is also an equivalent problem. All the results 
shown in the present paper may therefore be applied to satellite systems. The con- 
sideration of planets revolving in opposite directions involves new expressions for 
canonical variables presented in Section 3. We adapt the disturbing function of the 
Beauge & Michtchenko model to retrograde motions in Section 4. In Section 5, we 
express the Hamiltonian for retrograde resonances. An orbital fit of the HD 73526 
planetary system involving a 2: 1 retrograde MMR was found consistent with the 
observational data in Gayon & Bois (2008a). As a consequence, in Section 6, we 
apply the analytical model to the HD 73526 planetary system ruled in such condi- 
tions and compare our numerical and analytical results. Finally, in Section |7] our 
results obtained by the analytical expansion are visualized by surfaces of section. 
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2 General expression of the Hamiltonian 

We suppose a star of mass Mq and 2 planets of mass in\ and m2 such as 
Mq » m\,m2 and moving around their barycenter. We choose to express the 
Hamiltonian of the system in the heliocentric frame by using the following De- 
launay variables : 





Li 


= Piyfimi 




Gi 


= U^l-e^i 


Qi 


Hi 


- GjCos/, 



(1) 



with ai,eiJi,Qi,(x)i, and M, the geometrical orbital elements of the planet i, and 
lii — G y/Mo + Mi. G is the gravitational constant. /?,■ = Momi/{Mo + irii) is the 
reduced mass for the planet i. 

The Hamiltonian of the system (denoted F) consists of a Keplerian part {Fq) 
and a disturbing function {Fi) : 



F - Fn + F] with 



i=i (2) 

1 miMo ^ 

Fi = -Gmim2- + -——{xiX2 + y\y2 + Z1Z2) 

A Mq 



The leading part of the disturbing function, called the direct part, depends on the 
instantaneous distance between the two planets (A). Owing to the choice of the 
origin of the coordinate systems, the second part of Fi is called the indirect part of 
the disturbing function. F\ depends on the barycentric velocities (i:,, j,, z,) of each 
planet, expressed in cartesian coordinates. 

To study the 3-body problem, different sets of canonical variables can be used, 
depending on the simplifications carried out for reducing the number of degrees 
of freedom. For instance, Beauge & Michtchenko (2003) applied the so-called 
modified Delaunay variables to solve the coplanar, prograde and eccentric 3-body 
problem. In the present paper, we define a new set of canonical variables, espe- 
cially devoted to planets revolving in opposite directions. 



3 Set of canonical variables for planets revolving in opposite directions 

In this study, we choose to set a prograde direction of revolution for the inner 
planet (noted 1) and a retrograde one for the outer planet (noted 2). From the 
previous set of Delaunay variables ([T|, we therefore define the two following sets 
of canonical variables for planets 1 and 2 respectively : 

. -wi = -(i3i + wi) ii-Gi = ii(l - a/i - <^\) (3) 

-Qy Gi-Hi = Li Jl-e]{l -cos/i) 
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' A2 = -M2 + cl>2 -L2 - ~ Pi y/fiiAi 

. -0)2 = -{Q2 - (^2) G2-L2 = - ^2(1 - a/i - (4) 

Q2 G2+H2 = L2yjl -ejil +COS/2) 

The first set of canonical variables (O corresponds to the well-known modified 
Delaunay canonical variables while the second one (|4]l is developed for retrograde 
planetary motion. /I, is the mean longitude of the planet i. We note that the longi- 
tude of periastron is defined by: oj = Q + o), for a coplanar and prograde planetary 
motion (which is equivalent to consider 7 = 0°) and hy ch - Q - oj, for a coplanar 
and retrograde planetary motion (equivalent to 7 = 180°). Restricting the study 
to a coplanar but counter-revolving problem, the set of canonical variables of the 
system can be simplified as follows : 

^1 U 

-Wi L\ - Gi 

~(1)2 G2 ~ 7^2 



4 Expansion of the disturbing function 

The method used for the expansion of the disturbing function originates from 
Beauge & Michtchenko (2003). They expanded an expression for the Hamiltonian 
that is valid for high eccentricities of one or both planets. Their method provides 
a good assessment of the disturbing function when planets are close to the colli- 
sion point. Hence, the solution converges in all points of the phase space except 
for the singularities (corresponding to planetary collisions). Moreover, the rate of 
convergence does not depend on the values of eccentricities but rather on the order 
of magnitude of the distribution function itself. Besides, depending on the studied 
planetary system, even for very eccentric orbits, the rate of convergence of the 
disturbing function may be relatively fast. 

Classical methods such as Laplace 1799 or Kaula 1962 were expanded for 
coplanar and quasi-circular orbits and were useful for studying dynamics of the 
Solar System (or hierarchical systems in general). Nonetheless, due to the de- 
tection of eccentric exo-planetary systems, new and specific expansions of the 
disturbing function were needed. As a consequence, the Beauge & Michtchenko 
method seems much more suitable for such eccentric planetary orbits. To apply 
this method for counter-revolving configurations, we propose to expand again the 
disturbing function, by considering the prograde motion of the inner planet and 
the retrograde motion of the outer planet. 



4.1 Direct part of the disturbing function 



The direct part of the disturbing function that generally encounters problems of 
convergence may be expressed as a function of the heliocentric radial distances r,- 
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of both planets and the ij/ angle between both bodies as seen from the star : 

^^{r\ + rl- In ri cos iliT^''^ (6) 

We note that for counter-revolving configurations, the angle if/ between both bod- 
ies is defined b}0 : 

<A = /l+/2+^^^ (7) 

with fi the true anomaly of the planet Acj - oji-cbi and W|i,2) defined in ([3]) and 
(O. A concise expression of the previous equation is written in ^ by taking into 
account the ratio p — rxjri such as : 

^ = (1 +p2-2pcosi/^)-i/2 

The key of the Beauge & Michtchenko (2003) method lies in the expansion in 
power series in a new variable noted x and corresponding to a measurement of the 
proximity of the initial condition to the singularity mljA : 

N 

^ = (l+x)-i/2^^M" (9) 

11=0 

with X - - 2pcosi/r. r2lA has a singularity at x - -1. The determination of 
the coefficients b„ is performed by the way of a linear regression for x values 
greater than - 1 + 6, 6 being a positive parameter close to zero. A good precision 
of the direct part of the disturbing function may be reached for a good compromise 
between the 5 value and the choice of A'^ order in the series expansion. Contrary to 
classical methods involving Fourier series of the ^ variable or power series in p, 
not only the convergence rate of this method is improved but also the expansion 
of the disturbing function can be applied for eccentric two-planet systems. More 
details can be found in Beauge & Michtchenko (2003). 

From ^ and by using the explicit expression of x, we find : 

Changing from powers of cos ^ to multiples of ijj and by using the explicit expres- 
sion of the t/f angle, one obtains : 

N N-l I >2h+/ / >-2«-/-l 

^.y y2A,,„«2-'P p cos(//,+//2+/^^) (11) 

With Au, = (-1) 2,=„ bu, (,^2,-2„) i r-„ )r/. y = I 1 if / > 

"* For both prograde orbits, the ifi angle is defined by: ip = fi - + with aii = Qi + aii. 
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The direct part of the disturbing function may be expressed in terms of mean 
anomaly by using the Fourier expansion of the following functions (e.g. Hughes 
1981) : 



Q"cos(//) = K;'cos{mM) 

in=—oo 

(^) sin(//) = Yj Kh^^imM) 



(12) 



with the Hansen coefficient function of the eccentricity (Kaula 1962) : 

oo 

- e\'-"'^Yj^:U^../' (13) 

.5=0 

with J'j^'^^j ^^^2 Newcomb operators, w\ - max(0, m-l) and W2 - max(0, l-m). 
Let us recall that the Newcomb operators obey to simple recurrence relations 
(Brouwer & Clemence 1961; Murray & Dermott 1999). 
Substituting ( 1131 ) into ( 112b . we obtain : 



oo oo 



cos(//) ^YuYu B„,,j,„,e-'cos(mM) 

/=0 m = -oa 

(14) 

sin(Z/) =^ ^ S„,,j>, sin(mM) 

j=Q ni--oo 

The direct part of the disturbing function is therefore expressed as follows : 

oo oo W 2A' 

^ - J] YjJli ^«'-')/2 ^iUKm,n «' c{ e\ cos(mMi + nM2 + lAtb) 

j.k=Q ni,n=-(xi 1=0 i=0 

u ^^^^ 

s l"W[ , s 1-1^2 



4.2 Indirect part of the disturbing function 

From (|2]i, the indirect part of the disturbing function is given by the function Ti : 

mim2 

T\ = — — (xiX2 + yiy2 + Z1Z2) (16) 

Mq 

where i:,, y,, and z, are the barycentric velocities of the planet expressed in carte- 
sian coordinates. In the case of coplanar orbits, the expressions of x, and j, take 
into account the direction of motion of the planet /. We define x, for a prograde 
motion and a retrograde one as follow^! : 



dxi dxi dMj 
dt dM: dt 



^—^Hj for a progade motion 

«, for a retrogade motion 

oMi 



(17) 



^ A similar equation is obtained for the expression of y^. 
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Consequently, considering a coplanar problem and planets revolving in opposite 
directions, the indirect part is equal to : 



GOT1OT2 .W2 

1 1 = a ' 

02 



d 



d / Vi \ d 



dM\\a\j dM2\a2] dM\\a\j dM2\a2} 



(18) 



By using ( fT4l) and considering the prograde motion of planet 1 and the retrograde 
one of planet 2, one obtains : 







a\ 




Zi - 








X2 _ 


(fl2 


02 




- 


(fl2 


02 





00 00 
y=0 m=-oo 

CO CO 

y'=0 /?i=-oo 
00 00 

ZZ' 

^=0 /;=-oo 
00 00 



A'=0 n=-oo 

Substituting these results into ( fTSi t, we find : 

00 00 2A' 

^. — Z Z Z^v 



i,k=0 m,n=-oo i=Q 



X mnB\^\j^,„B\^\^k^ne\e\ cos(mMi + nM2 + zfw) 
where A, are constant coefficients such as : cv"'^' = Z^^^iQ"'- 



(19) 



(20) 



By combining ( 1151) and ( 1201) . we obtain the complete expression of the disturb- 
ing function as follows : 

Gin M °° 00 N IN 

F\ = ^'V,<:,«!,n,' o' e\ e\ cos(mMi + nM2 + lAcb) (21) 

J.k=0 m.n=-oa 1=0 i=0 

where = Aia-DiiDij^jXm,,, - <5/,i AimnBi,ui„,Bijxn are constant coeffi- 
cients. Rij,k.im,n.i are independent of initial conditions and then require to be de- 
termined once. We are now able to express the angles of resonance in the case of 
counter-revolving configurations as shown in the following section. 



5 The resonant average Hamiltonian 

Considering planets close to a MMR and revolving in opposite directions, we 
define the MMR ratio hy p + q/p with p + Q and p < 0. For instance, when two 
planets revolve in opposite direction and have a period ratio of 2, the q order of 
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resonance is equal to 3 and p - -1. We set s - p/q and define the following set 
of canonical variables in the case of a retrograde MMR : 



h J2^-L2-{\+s){h+h) 



o-i = (1 + s)A2 - sAi -0)1 /i = Li (1 - yj\-e\) 

(T2-{\ + s)A2 - sAi -0)2 h - - L2{\ - yjl - el) 



(22) 



where cr\ and o"2 are the resonant angles while /i and I2 are the conjugate mo- 
menta depending on the eccentricities e\ and 62. 

Let 6 be the following angle of the disturbing function such as : 9 - mM\ + 
nM2 + lAo). Considering the new set of variables, the expression of 6 becomes : 

- mcTi - n(T2 + /(a"2 - cri) + [in(p + q) - np]Q (23) 

where qQ - {Ai -A2) is the synodic angle. Hence, the disturbing function depends 
on three angular variables (cti , 1X2, Q). Considering the new set of canonical vari- 
ables (0-1,0-2, Q, A2, I\,l2,qJi,J\ + J2), we find the following constant of motion : 

Ji+ J2- constant (24) 

As a consequence, our system has therefore two integrals of motion: the F 
Hamiltonian and J,ot — J\ + /2- It is well known that the frequency of the angle 
Q is much higher than that of o-, . As a consequence, we consider only long period 
perturbations in the Hamiltonian. The system is then averaged with respect to the 
synodic angle as follows : 

-In 

FidQ (25) 



--\ 

27r Jo 



The averaging over the Q variable written in ( |25] l implies a third constant of mo- 
tion, namely 7] . As a consequence, the system consists of four degrees of freedom 
and three constants of motion. From now on, the system can be reduced to four in- 
dependent variables, namely cr\ , 0-2, /i and 12 - Taking into account the condition of 
MMR (« - m{p + q)lp), we obtain the final expression of the disturbing function 
I and the average Hamiltonian ( |27| l respectively, for a retrograde resonance : 



Gm]m2 



k,m,x m,„,- /„„,, IN 



'^2 j^Q m=-m,„„ 1=0 i=0 (26) 

X tt'CjCj cos((ot - /)o-i + (/ - n)a-2) 

/"^/^^ GmiOT2 Y"' - i k 

^ ^ ~ Zj T/T Zj ^''j-'''"'-' ^1 ^2 cos((m - /)o-i + (/ - n)o-2) (27) 

■ 1 J^I-'; ^2 ■ - / I 

Due to D'Alembert's properties of the disturbing function, some coefficients 
are null if one (or more) of the following conditions is reached : 1) j < \m- l\,2) 
k <\l - n\, 3) (m - I) even (odd) number and j odd (even) number, 4) (Z - n) even 
(odd) number and k odd (even) number. 
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6 Comparison with numerical methods 

Using the final expression of the Hamiltonian, the Ji and J2 constants of motion 
as well as the following set of canonical variables. 



(Ti - {I + s)A2 - sAi - a)i h - Liil-Jl-eh 

I (28) 

0-2 = (1 + s)A2 - sAi -a>2 4 = -^2(1 - a/i - ^2) 

we can firstly integrate a two-planet system both in counter-revolving configu- 
ration and in MMR and secondly, compare the analytical results with numerical 
ones. Using a Bulirsch-Stoer method, we integrate numerically the following dif- 
ferential equations : 



OF _ dFi 

d^i- (29) 
OF dF\ dFo 

For the integration of (|29l l, we use our initial conditions for the HD 73526 
system, located very close to the 2:1 retrograde MMR (i.e. at the edge of a V- 
shape structure in a stability map in [a, e] orbital elements; see Gayon & Bois, 
2008a) : 

Mo = 1.08 Mo 

Ml ^ 2.9 Mj,p M2^ 2.5 Mj,p 

fli=0.66AU 02 = 1.05 AU 

(30) 

ei=0.19 62=0.14 ^ ^ 

0-1 = 94 (deg) 0-2 = 94 (deg) 

Such initial conditions (130b are sufficiently close to the retrograde MMR to ap- 
ply our analytical expansion^ These initial conditions (more paticularly the semi- 
major axes) are averaged for the analytical expansion. 

Fig-fflshows the time evolution of eccentricity of both planets. Dots represent 
the numerical solution while black curves the analytical oneQ Both solutions ex- 
press the same behavior, with a relative error of 1.2% and 10%, on average, over 
the eccentricity of the inner orbit and the outer orbit, respectively (absolute er- 
rors being 0.002 for ei and 0.005 for 62 on average). Due to the averaging of the 
Hamiltonian over short periods, numerical dots are scattered on both sides of the 
analytical solution. 

Fig.|2]shows the time variation of the variable Acr. Because the 62 eccentricity 
periodically reaches the zero value, the numerical method does not always permit 
to determine the value of the Acr angle. This phenomenon is expressed by the 
vertical scattering of dots from to 360 degrees. As a consequence, surfaces of 

^ We forecast to study various initial conditions in a forthcoming paper, notably for a "full" 
retrograde resonance. 

' More details on the used numerical method and corresponding results may be found in 
Gayon & Bois (2008a). 
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section obtained numerically comprise this numerical bias. The analytical method 
is therefore more reliable when the eccentricity of a planet reaches values close to 
zero. Besides, the analytical approach ensures the properties of stability found for 
each planetary system. 



(a) 0.2 I < : 1 1 1 (b) 0.16 




300 600 900 1200 1500 300 600 900 1200 1500 

Time (years) Time (years) 



Fig. 1 Time variation of the eccentricites ei and e2 for the HD 73526 planetary system in 2;-l 
MMR. Red dots represent the numerical solution while black curves the analytical one. Used 
initial conditions are given by: Mo = 1.08 Mq ; Mi = 2.9 Mj„p ; M2 = 2.5My„p ; ai = 0.66 ; 
612 = 1.05 ; ei = 0.19 ; £2 = 0.14 ; cr, = 94 ; 0-2 = 94. 




300 600 900 1200 1500 
Time (years) 



Fig. 2 Time variation of the Act variable (where Act = cti - 0-2) for the HD 73526 planetary 
system in 2:- 1 MMR. Red dots represent the numerical solution while black curves the analytical 
one. Used initial conditions are the same as in Fig.[T] 
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7 Surfaces of section 

Nature of planetary systems can be directly deduced from the behavior of the 
resonance variables. As a consequence, from the previous Hamiltonian expansion 
( 127b of the three-body problem, we can plot surfaces of section and study the 
dynamics of two planets in retrograde MMR. 

Fig. [3] shows surfaces of section plotted for the HD 73526 planetary system 
previously defined. The energy level corresponding to the initial conditions i30i 
was numerically evaluated as -0.13835250 (units of AU, solar mass, and year). 
With our analytical expansion, we obtain the Hamiltonian value of -0.13835197, 
which is in good agreement with our numerical results. Panel (a) of Fig. [3] cor- 
responds to the inner planet and is plotted in the {e\ cos ct\ , e\ sin cri) parameter 
space. The section plane for the inner planet is chosen such as (T2 - 0. Similarly, 
panel (b) is plotted for the outer planet in the (e2 cos crj , si C2) parameter 
space and using the plane cti - 0. The initial conditions corresponding to ( [30l) are 
plotted in red. All the solutions found for F = -0.138352 are quasi-periodic. 

In panels (a) and (b) of Fig. [3] we represent the circulation state of the ci and 
0-2 variables in black and red colors. The libration state about 0° is plotted in blue; 
such initial conditions are located inside the 2 : 1 retrograde MMR. Nevertheless, 
for all the sets of initial conditions required for plotting Fig. |3] the time variation 
of o"! and cr2 expresses a circulation (in time) of both variables (not shown here). 
This time circulation seems contrary to the possible behavior in libration of ct] 
and cr2 shown on our surfaces of section. Such initial condition sets inside the 
2 : 1 retrograde MMR and their behaviors could therefore point out a particular 
characteristic of retrograde MMR. Such a characteristic would deserve a specific 
study. 





^9 







(b) 0.15 



-0.2-0.15-0.1-0.05 0.05 0.1 0.15 0.2 
e, cos o, 




-0.15 



-0.15 -0.1 -0.05 0.05 0.1 0.15 



Fig. 3 Surfaces of section for the inner (a) and the outer (b) planets of the HD 73526 system 
found close to a 2:-l MMR. The value of F is -0.138352. Initial conditions l l3Qt are represented 
by the red curve. Black and blue curves respectively correspond to circulation and libration of 
the IT, variables in the plane crj^i = 0. 
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8 Conclusion 

In the present paper, we investigated the three-body problem in the particular case 
of retrograde resonances. Our study is derived from the Hamiltonian approach of 
Beauge & Michtchenko (2003), which was expanded for eccentric, coplanar and 
prograde orbits. To apply this method to the retrograde MMR case, we expanded 
again the disturbing function, when considering a prograde motion of the inner 
planet and a retrograde motion of the outer planet. Hence, we defined a new set of 
canonical variables, which allow us to express correctly the angles of resonance 
in the case of counter-revolving configurations. 

Although the exploration of the A^-body problem is accessible from numerical 
methods, the acquiring of an analytical "rail" notably contributes to a deeper un- 
derstanding of the numerical investigation. As shown in Section |6l the analytical 
method also permits a better determination of the resonance variables. Moreover, 
from the behavior of the resonance angles displayed in surfaces of section, we can 
directly infer the local behavior of a 3-body system, that is to say its stability or 
its chaoticity. 

Until now, no planetary system has been truly detected in counter-revolving 
configuration. However, given the efficiency for stability of retrograde MMR (see 
Gayon & Bois 2008b), such a detection might occur. The work presented in this 
paper is firstly based on a dynamical study of theoretical two-planet systems. Nev- 
ertheless, a dynamical system composed of two satellites orbiting a planet is an 
equivalent problem. Since some satellites of Saturn are found in counter-revolving 
motions, our expansion of the three-body problem solved in the case of retrograde 
motions could then be applied for satellites of the Solar System. 

Acknowledgements We thank the anonymous referees for their useful comments and sugges- 
tions that helped to improve the paper. 



A Appendix 

This appendix shows the analytical expansion when the body moving on a retrograde orbit is 
the inner planet. The new set of canonical variables for the resonant averaged Hamiltonian is 
written as follows (equivalent tol22b : 

M Ji=-U+ s(Ii + h) 

A2 J2= L2-{\ + s){h+l2) 



o"! = (1 + s)A2 - sAi - (h\ /i = - Li (1 - -yjl-ej) 

0-2 = (1 + s)A2 - sAi -0)2 h = L2 (1 - -^1 - ej) 



(31) 



with Ai = -Ml + Cj\, A2 = -M2 - c1j2, ciii = Oi - oji and 012=02 + (^2- 

Although the set of canonical variables changes, the expression of the resonant averaged 
Hamiltonian remains the same : 



2 ,,2o3 



Zfijfij Gmim2 v - it 
2j Ri.j.kjn.iaie\e2C0!i((m-l)cri +(1-11)0-2) (32) 



;=1 ( tj.k.iuj 
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